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Abstract 

Based on the BCS theory we outhne a conventional path integral derivation 
of the transverse force and the friction for a vortex in a superconductor. 
The derivation is valid in both clean and dirty limits at both zero and finite 
temperatures. The transverse force is found to be precisely as what has been 
obtained by Ao and Thouless using the Berry's phase method. The friction 
is essentially the same as the Bardeen and Stephen's result. Q 
1. Introduction 

In condensed matter there are two types of important excitations: the gapless and el- 
ementary ones, such as spin waves or phonons, which determine the local properties, and 
the topological ones, domain walls or dislocations, which determine the global properties. 
Vortices in superconductors belong to the second category. They have been under intensive 
studies since earlier sixties0. Still, there has no general agreement on form of the vortex dy- 
namics equations. We point out here that the main disagreement among many microscopic 
derivations lies in the use or not of the relaxation time approximation to account for the 
effects of impurities in the force-force correlation function. To provide a natural treatment 
of the both the transverse force and the friction and to avoid the use of this approximation. 



^ Supported by Swedish NFR. To appear in Physica C. 
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we present a path integral derivation within the BCS theory. The important results here are 
that the transverse force is the same as obtained by the Berry's phase method! and by the 
total force-force correlation function methodi, and is indeed insensitive to impurities, and 
that the friction coincides essentially with the Bardeen and Stephen result^ but differs in 
some important details. Specifically, the nature of the friction is determined by the spectral 
function of the Hamiltonian. The relevance of the present theory to experiments has been 
discussed elsewhere .@ 
2. Path Integral Formulation 

We begin with the standard BCS Lagrangian in the imaginary time representation for 
s-wave pairing. We will only consider neutral superconductors here, but the coupling to 
electromagnetic fields does no affect our results in the extreme type II case. The Lagrangian 



IS 



Lbcs = EcT i^l (a;, r) ( - jJp 

-i^v' + vix))M^,T) (1) 

-9^\{x,r)tljl{x,T)^l{x,T)^i;{x,T) , 

where ipa describes electrons with spin a = (t, i), fip the chemical potential determined by 
the electron number density, V{x) the impurity potential, and x = {x, y,z). A vortex at x^ 
has been implicitly assumed. The partition function is 

z= /DK,v^V'}x 

exp{-y,^^dTjd'xLBcs} , 
with P = l/ksT , and d'^x = dxdydz. Inserting the identity in the functional space, 

1= JV{A*,A}exp{-^J,^^dTjd'xx 
^^^1+ iA(x,r) 

into Eq.(2) we have 
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Here the Hamiltonian is defined as 



H(A, A* 



V 



(3) 



H A 

with if = —{h^/2m)V'^ — fip + V{x). Integrating out the electron fields ip^ and first, 
then integrated out the the auxiliary(pair) fields A under the meanfield approximation, one 
obtains the partition function for the vortex 



with the effective action 



Z = JV{x4exp{-^} 



Tr In - ^ Jq'^ dr J d^x\A\ 



(4) 



(5) 



where Tr includes internal and space-time indices, and the Nambu-Gor'kov (NG) Green's 
function G defined by 



{hdr + n)G{x, t; x\ t') = 5{t - t')5\x - x'). 



together with the BCS gap equation, or the self-consistent equation. 



(6) 



A(x,r) = -g < V'|(x,r)^-f(x,r) > . (7) 

A special attention should be paid to the equal time limit of the NG Green's function.0 

We assume that the vortex is confined to move in a small regime around a point at xq, 
which allows a small parameter expansion in terms of the difference between the vortex 
position Xy and xq. We look for the long time behavior of vortex dynamics under this small 
parameter expansion. For the meanfield value of the order parameter, this expansion is 



(8) 



Here 5xy — x^ — xq. In Eq.(8) we have used the fact that when x^ — xq the vortex is static. 
The effective action for the vortex to the same order is, after dropping a constant term, 



with 



5x^ ■ VqAq 5x^ ■ VqAo , 



^ 8xv ■ VqAo ^ 



(9) 



(10) 



\^ 5x^ ■ VoA* 

Here Go is the NG Green's function with A(A*) replaced by Ao(Aq), and the gradient Vo 
is with respect to Xq. 

Now we construct the NG Green's function Go following the usual procedure. First, we 
consider the eigenfunctions of the Hamiltonian Tic = H(Ao, Aq). The stationary Schrodinger 
equation, the Bogoliubov-de Gennes equation, is 



(11) 



with 



( ( A 



Since Hq is hermitian, all its eigenstates form a complete and orthonormal set, that is. 



Jd^x'^i{x)^k'{x) = 6k,k' , 



and 



6^{x — x') , 



with '^\x) — {u*{x),v*{x)). Furthermore, Eq.(ll) has an interesting property: 



with 



^ v*{x) ^ 



V 



-u*{x) 



Finally, since Sx^ ■ VqHq — E', for k ^ k' we have 

! d?x-^l{x)T.'-^k'{x) 

= - Ek)5x^ ■ jd!^x'^l{x)Vo^k'{x) , 

which implies that the two ways of specifying the vortex coordinate, through the trapping 
potential or through the order parameter, are equivalent. Those relationships will be used 
below. 

Given the eigenfunctions of Eq.(ll), the NG Green's Go can be expressed as 

Go{x, r; x', r') = En,, g'SS?^^!^)^!!^')- 

Here cUn = mr/hp, with n odd integers. Assuming the rotational symmetry after the impurity 
average, a straightforward calculation leads to the following effective action 



Seff= ljfdT\K\Sx,{T)\'' 



+ jf dT'F\\{T - t')\5x,{t) - 5x,{t')\^ 
+ dr'F^ir - t'){5x,{t) x 5x,{t')) ■ z 

with the spring constant in the effective potential. 



(12) 



(13) 



the longitudinal correlation function. 



(14) 



and the transverse correlation function, in the long time limit, in terms of the virtual tran- 
sitions. 



F±(r) = -dr-r'S{T-T')Ek,k'Id^xJd^x' 

Hfk - fk')\z ■ [ml{x')Wo<i>u'{x')x (15) 

or in terms of the contribution from each state, 

F±(r)= -dr-r'S{T-r')Y.kSd!'xz- 

hifkVoulix) X VoUkix) (16) 
-(1 - fk)Vovl{x) X VoVkix)) . 
To reach Eqs. (12-16), following two identities have also been used: 



g — s"n<5 




Q- 



and 



_1 e-»(^n-u.„/)T ^ l cosh[f (M_|^|)] 
^n-"' /3 ih{ur,-ic„,)-E 2 sinhfE^l 



2 J 



with the Fermi distribution function = 1/(1 + e^'^*=), and the spectral function 

J{^) = jEk,k'Sihuj - \Ek- Ek'Dlfk- fk'\x 



(17) 



\jd^x¥,{x)Vono^k'ix)\' 

In the following we discuss the implications of K, F\\, and Fj_ one by one, and show that 
Eq.(12) contains both the dissipative effect corresponding to the Bardeen-Stephen resultB 
and the transverse force identical to the one by Berry's phase methodll or the total force 
correlation function method!. 
3. Correlation Functions 

3. 1 Effective Spring Constant K 

For the purpose of getting the friction and the transverse force, the precise value of K is 
irrelavent. 

3.2 Longitudinal Correlation Function F\\ 

This longitudinal correlation function contains all information on the vortex friction, 
revealed by the fact that Eq's.(12,14) are identical to the time-nonlocal action in the quantum 
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dissipative dynamic^. The friction is determined by the low frequency behavior of the 
spectral function J{uj), Eq.(17). The temperature is naturally incorporated into the spectral 
function. This form of the action has been used in the study of the vortex tunneling.0 

If we classify the eigenstates of the Bogoliubov-de Gennes equation, Eq.(ll), according 
to core (localized) and extended states, the dissipation comes from all three parts: core- 
to-core, core-to-extended (or extended-to-core), and extended-to-extended contributions. 
Those contributions are added up, because they are all positive, which differs from the 
situation for the transverse correlation function. 

For a clean and neutral superconductor, the core states are separated by finite energy 
gaps from each other. Their contributions to the dissipation will be exponentially small 
in the slow motion case. The dominant contribution in this case comes from the extended 
states, because their energy spectrum is continuous. This is also true for a thin film after 
the consideration of the electro-magnetic field effects. Those contributions correspond to 
super-Ohmic cases.li^ 

The presence of impurities mixes up all the clean limit eigenstates, and generates a 
continuous distribution of the energy, even for the core states, after the impurity average. A 
perturbative calculation shows that it is indeed the case, and the spectral function becomes 
Ohmic.0 This suggests that in the clean limit the dissipation is super-Ohmic, and turns into 
Ohmic in the dirty limit. Hence in the dirty limit one should expect the Bardeen-Stephen 
result!, because they have treated the core as in a normal state, corresponding to the uniform 
distribution of the energy eigenvalue. 
3. 3 Transverse Correlation Function F± 

In the long time limit the there are two equivalent forms for the correlation function: 
The virtual transition expression, Eq.(15), and the individual state contribution, Eq.(16). 
Their equivalence has already been discussed in Ref.i with the aid of conservation laws. 
Eq.(15) can be again casted into various equivalent forms because of the cancelations among 
those virtual transitions. For example, it can be expressed only in terms of the core-to-core 
transitions. In this case it is a statement of the spectral flowi. But one can also cast it into 



the forms of core-to-extended, or extended-to-extended transitions. 

In Eq.(16), the counting of individual state contributions is expressed as an area integral 
of the momentum commutator. If we express it as the line integral of the momentum density 
far away from the core, the immediate conclusion, as drawn in Ref.il, is that localized core 
states do not contribute, because the momentum density at the core is zero. The insensitive 
to the impurities is most transparent from Eq.(16): In the one-body density matrix both 
the electron number density and the phase 6{x all insensitive to impurities. Here 

the phase 9 is defined through the order parameter A(x, — \A\e'''^^^^~^'"\ with q = ±1 
describing the vorticity along the z direction and 6{x) = arctan(?//x). Hence, as long as 
the localization effects due to impurities are negligible, such as in the cases of usual dirty 
superconductors, the transverse force will not be influenced by impurities. An analogy can 
be drawn to electrons in semiconductors and in metallic rings: The momentum is insensitive 
to band structures and disorder potentials. 
4. Discussions 

1. The major difference between the present approach and those in Ref 's.i0 lies in the use 
or not of the relaxation time approximation in the calculation of the force-force correlation 
function. It is found in Ref's.i'0 that this approximation is needed to makes the friction 
finite, which at the same time reduces the magnitude of the transverse force. Early studies 
of this approximation have shown that it is in general ambiguous and problematic in the 
force-force correlation functional. The source for the ambiguity is that it is usually unable 
to know a priori the transverse, the friction, and the fluctuating forces, although the total 
force is often well defined. This is exactly what happens in the present situation. 

In the present approach without this approximation we are able to obtain the friction 
determined by the spectral function of the Hamiltonian, and demonstrate that the transverse 
force is insensitive to impurities. 

2. The explicit counting of each state contribution, Eq.(16), confirms that the contribute 
to the transverse force comes from a region away from the core. One way to show it, as 
discussed above, is from the consideration of the momentum density, which is zero at the 
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phase singular point, the vortex position. In an equivalent way, if the one-density density 
matrix, defined by {uk, vl}, is finite at the vortex position, as allowed by the Bogoliubov-de 
Gennes equation, the corresponding Uk or will not have the phase factor 0.0 Hence, there 
is no an additional contribution to the transverse force from the phase singular point as 
supposed in Ref.i. 
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